The role of controllability in optimizing quantum dynamics 



Re-Bing Wu 

Department of Automation, Tsinghua University, Beijing, 100084, China 
Center for Quantum Information Science and 
Technology, TNList, Beijing, 100084, Chinch 

Michael A. Hsieh 

Department of Chemistry and Center for Quantum Information Science and Technology, 
University of Southern California, Los Angeles, CA 90025, USA 

Herschel Rabitql 
(Dated: January 7, 2011) 

Abstract 

This paper reveals an important role that controllability plays in reducing the complexity 
of optimizing quantum control dynamics. Based on the quantum control landscape, which 
refers to the optimization objective as a functional of the control fields, we analyze the critical 
topology of the gate fidelity landscape for uncontrollable systems due to SU{2) and SU(3) 
dynamical symmetries. The loss of controllability leads to multiple local suboptima on the 
control landscape, which form traps on the way towards seeking the global optimal value for the 
objective. The control landscape generally becomes rugged when controllability is constrained 
by reduced dynamical symmetry, and even more rugged if in addition the target gate is not 
realizable. These results imply that the optimization of quantum dynamics can be seriously 
impeded by local search algorithms under these conditions. 
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I. INTRODUCTION 



A quantum system is referred to as controllable if it can be transferred from an ar- 
bitrary state to another through engineering the dynamics. As a fundamental aspect of 

'or various types of sys- 



1 11 ] closed systems and 



quantum control theory, controllability has been widely studied 
terns including infinite-dimensional (5), and finite-dimensional 
Markovian open systems However, the practical necessity for controllability is still 
unclear as it may not be required in many realistic applications and may not always be 
available. For example, the goal of realizing only a specific unitary transformation related 
to achieving a quantum computation algorithm, or reaching particular excited states of 
a molecule generally do not require full controllability. No apparent advantage can be 
seen by merely gaining additional controllability. The study of this issue is also motivated 
by the growing number of successful adaptive control experiments utilizing shaped laser 
pulses as the controls [2]. In most of these systems, there is no assessment of controllability 
and a simple yes-or-no controllability answer, likely has little practical relevance. 

Here we study this matter in the framework of the quantum control landscape 9[, which 
was developed for understanding the complexity of optimizing quantum dynamics [2J, [8J. 
A quantum control landscape is formally defined as a target functional, </[{£&(■)}], to 
be maximized (minimized) over the set of implementable control fields {ei(t), ■ ■ ■ , e m (t)} 
that steer the following iV-level quantum control system: 

dU{t) 



dt 



ih)- 1 \H + Y^t k (t)H k ^U(t), 17(0) = I, (1) 

where U(t) e U(N) is the system propagator at time t, Ho is the internal Hamiltonian 
and {Hi, ■ ■ ■ , H m } are the control Hamiltonians. 

A major goal of quantum control landscape theory is performance of a topological 
analysis of the critical points (i.e., the critical topology) from which insights can be gained 
about the complexity of optimizing quantum system dynamics. Consider the class of 
control landscapes in the form of J = J(U(tf)), where the system propagator U(tf) at 
t — tf is an implicit functional of the control field through the Schrodinger equation ([1]). 
The control landscape can be displayed either on the space of the control fields (i.e., in 
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the dynamical picture) or on the set of realizable propagators U(tf) (i.e., in the kinematic 



picture). In prior studies 
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131 ]. we have shown that the main topological characteristics 



of the critical points in the dynamical picture are preserved when mapped onto the critical 
ooints in the kinematic picture except for what appear to uncommon singularity issues 



17| . Thus, almost all maxima, minima and saddle points in the dynamical picture should 
correspond to the same type of critical points in the easier to analyze kinematic picture 
without distortion of the topological features. 

The kinematic control landscape may be constructed from the (low-dimensional) set of 
realizable propagators. If the system is controllable, then the set of realizable propagators 
form the unitary group. This circumstance has been studied for practical optimization 
objectives (e.g., the state transfer probability 9j, ensemble average 14| and gate fidelity 



10|). which collectively show that none of these landscapes possess local suboptima that 



may impede the search for the ultimate global optima. This behavior provides a basis to 
explain the large and growing body of experimental quantum control successes, although 
in the laboratory many ancillary issues and constraints may be present. 

Realistic quantum systems may not be always fully manipulable. A typical circum- 
stance is when the system is under some dynamic symmetry that restricts the system's 
evolution. Even if the system is not subject to a dynamic symmetry, the restrictions on 
available control resources (e.g., the amplitude, bandwidth and center frequency of the 
applied field) can still cause a practical loss of full controllability, in which case movement 
over the landscape when searching for an optimal control can become limited on the oth- 
erwise trap free landscape. Whether these uncontrollable systems still have a physically 
attractive landscape topology becomes practically important for assessing the complexity 
of searching for optimal controls, which is explored in this paper. 

The loss of controllability poses restrictions on the set of reachable states, which im- 
pacts the kinematic analysis. The set of reachable states for control-restricted systems is 
generally very complex to identify, while that for systems with dynamic symmetry may 
be simply revealed (e.g., as a Lie subgroup of U(N)). The landscape topology for systems 
with dynamic symmetry has been identified for gate fidelity under certain special symme- 



tries 



% ll5| , where no false traps were observed. However, as will be seen in this paper, 
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false traps (i.e., yet appearing as real in practice) can occur in the landscape of symmetry 
induced uncontrollable quantum systems. The paper will be organized as follows. Section 
II discusses the general conditions for the existence of critical points on such landscapes. 
Section III analyzes the case of the fidelity control landscape of transformations with 
SU{n) dynamic symmetry, where for n = 2 or n = 3 we derive the analytic structure 
for simple cases and graphically compare the landscapes for different symmetry types. 
Finally, conclusions are drawn in Section IV. 

II. GENERAL CONDITIONS FOR CRITICAL POINTS ON A LANDSCAPE 
WITH DYNAMIC SYMMETRY 

Let g be the Lie algebra spanned by (ift) _1 ifo, (ifi)~ 1 Hi, . . . , (ift) -1 if m , then the quan- 
tum propagator must evolve within the Lie group Q = exp g generated by g, which forms 
the dynamic symmetry of the system. If g forms a proper Lie subalgebra of u(N), the 
system must be uncontrollable. However, owing to the compactness of the unitary group, 
the set of realizable propagators can fill up the group Q if the control is not limited and the 
evolution time is sufficiently long. Hence, in a kinematic analysis, one can suitably study 
the landscape of symmetry induced uncontrollable quantum systems on the corresponding 
symmetry group. 

The condition for U e Q to be a landscape critical point can be derived through pa- 
rameterizing its neighborhood in Q as Ue A , where A e g represents the local coordinates. 
The variation of the landscape function with respect to A can then be written via the 
chain rule 

SJ g = (VJ(U), U5AJ = (u^J(U), 5A^ = 0, V 8 A e g, (2) 

where V J(U) is the gradient function over U(N) and the matrix inner product (X, Y) = 
ReTr(X^Y). Let B(U) = WVJ{U) G u(N), then the condition for U e Q to be critical 
is that 

(D(t/),g) = 0, (3) 

which requires that the vector D(£7) e u(N) be orthogonal to the subspace g. For example, 
D(£7) = [UpW,9] corresponds to the ensemble control landscape J = Tr(UpW6) and 
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D(Z7) = W'U — WW corresponds to the gate fidelity control landscape J = ReTr(W'U). 
The set of critical points with Q dynamic symmetry can then be expressed as 

M g = {Ueg-. (O(L0,g) = 0.} 

The condition ([3]) provides a clear geometric picture for the critical points, however, it 
is difficult to use for directly determining the critical topology with an arbitrary dynamic 
symmetry, except in the extreme case that the system is controllable (i.e., g = u(N)), 
where the condition can be simplified to D({7) = on U(N). It is also possible, as will 
be shown below, to simplify the analysis by making use of the symmetric form of the 
landscape function, from which we can either solve for the critical topology or easily 
visualize it for low dimensional symmetries. 



III. QUANTUM GATE FIDELITY CONTROL LANDSCAPE WITH SU{n) DY- 
NAMIC SYMMETRY 

The quantum gate fidelity landscape is defined as 

J(U) = N- l \Ti(W^U)\ (4) 

where N is the dimension of the system and W e U(N) is the target unitary transfor- 
mation [ljj. If the target transformation W is realizable, i.e., W G Q, then we can use 
the transformation U = W^U to change the landscape function into a canonical form 
Jg(U) = N- l \Ti(U)\. 

From group representation theory, a dynamic symmetry group Q C U(N) can be taken 
as a N- dimensional unitary representation of the abstract Lie group Q, which thereby 
can be decomposed into smaller unitary irreducible representations (UIR), i.e., one can 
always find a constant unitary similarity transformation in Q that turns each U G Q into 
a block diagonal form: 

( UW \ 

{ uv> J 
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where G U(Nj), Ni~\ — • + N r = N, and each sub-block cannot be decomposed 
into smaller ones. For the gate fidelity control landscape, the landscape function with Q 
carrying a reducible representation can be written as the sum of the landscape function 
values evaluated on its irreducible components. Hence, it is convenient to start with the 
case that Q is a unitary irreducible representation of some abstract Lie group in U(N). 

A classification scheme of uncontrollable systems was proposed 7J based on the Lie 
algebra g and their unitary irreducible representations (UIR). Here we only consider the 
class of SU(n) dynamic symmetry (i.e., the group of nxn unit trace unitary matrices). For 
example, the SU{2) dynamic symmetry may arise when the free and control Hamiltonians 
of spin-j particles (or a rigid rotor whose angular momentum is j 2 ) in a magnetic field 
are a linear combination of the spin operators J x , J y and J z that possess an N = 2j + 1 
dimensional UIR (denoted by T>j). Similarly, two non-interacting spin particles with spin- 
ji and spin-j 2 in a magnetic field (e.g., two electrons in a multi-electron molecular system) 
possess a tensor product of SU (2) UIRs V = <g> Vj 2 , which can be decomposed into 
the direct sum of smaller UIRs of SU{2). 

Suppose that Q carries an N- dimensional UIR of SU{n), then there must be a home- 
omorphic mapping 1Z that transforms any unitary matrix U G Q C U(N) to an n x n 
unitary matrix U' G SU(n) such that U = 1Z(U'). Let the complex numbers ei, • • • ,e n be 
the eigenvalues of U' corresponding to U, which satisfy |ei| = ■ ■ ■ = |e n | = 1 = e± ■ ■ • e n , 
then according to the well-known Weyl formula jy], the UIRs of SU{n) are labeled by a 
group of integers r\ > ■ ■ ■ > r n _i [6|, and the trace of U can be expressed as the following 
character function 



T±(U) = Xfi,-,r„_i(ei,-- • ,e r 



,e,l\ 



where 
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This formula shows that the original landscape function on an (iV 2 — l)-dimensional group 
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manifold Q can be reduced to a function of n— 1 independent variables (e 1; • • • , e„_i) on the 
(n— l)-tori T n_1 , thereby all U G corresponding to a critical point (ej, • • ■ , e n -i) £ T™^ 1 
form a subset of the critical submanifold with the same J value. 

As a direct application of the above formula, the case of SU(2) can be computed. 
Suppose that Q carries a UIR Vj of SU(2), then the landscape function is 



2j+l _ 2j+l 

J Vj (e u e 2 ) = (2j + l)- 1 ^ 2 — (5) 

ei - ^2 

according to ([5]). Let ei = e l/3 and e2 = e~ l/S , where (3 G [0,27r), then (j^J) can be reduced 
to a single variable function 

_ sin [(2j + 

J ^ (/3) - (2 J + l)sin/3 ' (6) 

whose critical points /3 2 , . . . characterize the critical topology. 

Besides SU{2), the SU(3) case is relatively simple, but a universal solution is still not 
available for the critical topology of all its UIRs. However, as the corresponding character 
formula involves only two independent variables, it is possible to solve for the critical 
topology of low dimensional UIRs, and, numerically, it is easy to visualize the landscape 
in a 3D image for any individual UIR (labeled by two integers ri > r 2 ). The dimension 
of the (r 1; r 2 ) UIR is given by N ri>r2 = DJ^ZlzlA_ T n ^he following, we will use SU{2) and 
SU(3) examples to study the impact of the symmetry group, the UIR dimensionality and 
the influence of realizability of the target gate upon the landscape topology. The case 
SU (n) for n > 3 will not be discussed in this paper. 



A. The dimensionality of UIR 

Consider the SU{2) case with W being realizable. Fig. 1 shows two examples with 
j = 3 and j = 7/2 for G [0, f] (the remaining part G [| ,2tt] is repeating and thus 
omitted) according to the formula respectively, corresponding to dimensions N = 7 
and N = 8. Both landscapes possess four local maxima, and their local minima are all of 
value zero. No saddle points are found. There are three traps on the landscape, and the 
results obtained from the two examples can be generalized as below: 



7 



Theorem 1 Denote \J\ as the largest integer that is no greater than j. The control 
landscape Jx> 3 has [j\ local suboptima. 

The systems that carry higher- dimensional UIRs (i.e., with a larger j in T>j) are less 
controllable because the number of degrees of freedom is always 3 (versus the total number 
(2j + 1) 2 in U(2j + 1)). Thus, we can see from the Theorem that the search for an optimal 
control is more likely to be trapped on such landscapes because there tends to be a larger 
of number of traps. Moreover, the convergence region of the global optimum U = I can 
be identified as f3 G [0, 2j+i\ according to (jSJ), which also shrinks when j increases, leaving 
the system more likely to be trapped by some of the — 1 local suboptima outside the 
region. Therefore, the loss of controllability leads to an increase of the ruggedness of the 
landscape. 




0.1 0.2 0.3 0.4 0.5 

P/Jt 

Fig 1: The character functions for P3 and T> 7 / 2 UIRs. 
B. The realizability of the target transformation 

Consider again the dynamic symmetry SU(2). If the target gate is not realizable, i.e., 
W ^ Q, then the dynamic symmetry group Q is not invariant under the transformation 
U = W^U, and hence the approach used above for realizable W is no longer applicable. 
In this case, the landscape function cannot be reduced to having a single characterizing 
parameter, and we have to analyze it on the original three-dimensional group manifold. 
Here we adopt the well-know Euler decomposition 
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which parameterizes the landscape function by (9,ip\,ip2). To facilitate visualization of 
the results, we consider the special case that the target gate W commutates with L z , i.e., 
[W, L z ] = 0, for which the landscape function reduces to the following two- variable form: 

J(U) = N- 1 \Tt{W^e ,hI "e ei "e^ I "}\ 

= jV-^Tr^e^+^^e^}! = J(6,<j>), 

where <p = ipi + ip 2 - F° r example, suppose that L z is diagonal under some chosen UIR 
basis, then the landscape function can be reduced to a two-variable function if the target 
gate is a diagonal iV-dimensional unitary matrix, whose diagonal elements affect the 
landscape topology. For such a class of unrealizable target gates, one can easily depict 
the two-dimensional landscape topology as exemplified in Figj2j Compared to the case of 
a realizable W, the landscape hight (i.e., the maximal fidelity) decreases to J max = 0.8 due 
to the inability of perfectly realizing the gate W. In addition, the landscape with realizable 
W possesses three critical submanifolds (owing to the intrinsic dependence upon a single 
variable), while the landscape with an unreachable W is much more rugged, implying 
that having an unattainable target state not only degrades the maximal achievable fidelity 
but also increases the complexity of the landscape topology. Consequently, local search 
algorithms (e.g., gradient based methods) on such a control landscape likely will be less 
effective. 

C. Comparison between different dynamic symmetry groups 

To compare the influence of different dynamic symmetry groups, we choose the (5, 2) 
and (6, 1) UIRs of SU(3) and the UIR V 7 of SU(2), whose dimensions are all fifteen (i.e., 
corresponding to different dynamic symmetries in the same physical system) . We directly 
show the 3D landscape profiles in Fig. [3j The landscape topology is less rugged for the 
<SW(3) all critical points are isolated. Moreover, the region of attraction for the 

global maximum is rather large, while that of the SU{2) symmetry is small, making its 
peak very sharp. This shows from another perspective that the enhanced controllability 
(corresponding to a larger dynamic symmetry group) helps to improve the search efficiency 
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Fig 2: Two-dimensional illustration of the landscape topology of J(U) = N 1 \Tr(W^U)\ with 
(a) W = I 8 € SU{2)- (b) W = diag{-l, 1, • • • , 1} £ 5W(2). 

for global optimal controls by smoothing the landscape topology. It is also interesting that 
the landscapes of different UIRs of the SU(3) symmetry group can be distinct (i.e., the 
UIR (5, 2) landscape is less rugged) even when their dimensions are identical. 

IV. DISCUSSION 

The present study on the landscape of uncontrollable systems with SU(2) and SU(3) 
dynamic symmetries show that controllability plays an important role in maximizing uni- 
tary gate fidelity, whether or not the desired target is attainable. The loss of controllability 
can lead to false traps over the corresponding restricted kinematic control landscape on 
the dynamic symmetry group. The ruggedness generally increases when controllability is 
weakened, i.e., either with a smaller dynamic symmetry group or an unattainable target 
unitary transformation, making the search for global optimal controls much more difficult. 

The details of the landscape topology depend on the particular structure of the dynamic 
symmetry group and the unitary representation it carries (see Q] for a classification of 
such systems). Uncontrollable systems with dynamic symmetry generally possess rugged 
control landscapes, but this may not always be the case. For example, the gate fidelity 
landscape with W realizable and the following reachable sets in U(N) were found (or can 
be easily proved) to be devoid of local suboptima: 
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Fig 3: The gate fidelity landscape topology for 15-dimensional systems under SU{2) and SU(3) 
dynamic symmetries: (a) the landscape with SU(2) symmetry with UIR P7; (b) the landscape 
with SU(3) symmetry with UIR labeled by (6, 1); (c) the landscape with SIA{2>) symmetry 
with UIR labeled by (5,2). Their landscape ruggedness is in decreasing order. 

1. Symmetry matrices in U(N) 4|; 

2. Self-dual matrices in U(N) |4|; 



3. Symplectic matrices (symplectic group) [151 ]; 



4. Real orthogonal matrices (orthogonal group) in IA(N). 

The landscape structure for the observable control problem [l(| also becomes rugged when 
the system is under dynamic symmetry. The corresponding the critical topology structure 
is dependent on the spectra of the system density matrix and the observable, and will be 
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more complicated than the case of unitary fidelity landscape. This important problem 
will be explored in future studies. 
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